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Abstract. We study spectra of alloy- type random Schrodingcr operators on 
metric graphs. For finite edge subsets of general graphs we prove a Wegner 
estimate which is linear in the volume (i.e. the number of edges) and the length 
of the considered energy interval. The single site potential of the alloy-type 
model needs to have fixed sign, but the considered metric graph does not need 
to have a periodic structure. The second result we obtain is an exhaustion 
construction of the integrated density of states for crgodic random Schrodinger 
operators on metric graphs with a Z'^-structure. For certain models the two 
above results together imply the Lipschitz continuity of the integrated density 
of states. 



In the present paper we study spectral properties of random Schrodinger opera- 
tors on a metric graph. More precisely, wc consider a countable metric graph and 
a random Hamiltonian with a so-called alloy type potential. Under suitable as- 
sumptions wc arc able to prove a Wcgncr estimate for the restriction of the random 
Hamiltonian to a finite part of the graph. Our Wegner estimate is optimal in the 
sense that it is linear in the energy and the volume. 

In the case that the integrated density of states (IDS) of the random Schrodinger 
operator exists, the Wegner estimate implies that the IDS is Lipschitz continuous. 
For certain random operators on Z'^-metric graphs we establish the existence of the 
IDS via an exhaustion procedure. More generally, the existence of an selfaverag- 
ing IDS is ensured by certain amenability assumptions on the metric graph and 
ergodicity assumptions on the random Hamiltonian, see Remark 4. 

The literature on Wegner estimates for Schrodinger operators in the continuum, 
i.e. on M"^, is quite extensive. Wc refer to the textbook accounts [CFKS87, CL90, 
PF92, StoOl, Ves04, KM] and the references therein. For the application to spectral 
localization in random media, a quite weak form of Wegner estimate is sufficient. 
However, upper bounds which arc linear both in the volume and the energy interval 
length are of independent interest, since they imply the regularity of the IDS. To 
obtain such a Wegner estimate, the proofs which are presented in the literature need 
some assumptions, beyond those necessary for a weak form of Wegner estimate. In 
particular, mostly the covering condition 



was assumed. Here u is a compactly supported, non-negative single site potential. 
Papers which are devoted to the question how this assumption can be removed or at 
least weakened include [Klo95, Kir96, Vcs96, CHNOl, KV02, Gic, CHK0:3]. In the 
recent [CHK] a linear Wegner estimate is proposed without the abovementioned 
covering condition. However, the deterministic part of the Schrodinger operator 
needs to be translation invariant. The reason for this assumption is that it implies a 
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uniform, translation invariant kind of unique continuation principle for the solutions 
of the Schrodinger equation, cf. Section 4 in [CHK03]. 

In our setting, on metric graphs, the situation is somewhat simpler, since in 
certain aspects the Schrodinger operators behave as in the one dimensional case. 
In particular, on the (one dimensional) edges one has a stronger (and simpler) 
unique continuation principle at disposal. We use techniques developed for Wegner 
estimates for operators on L^(K) in [KV02]. (See also [V<>s96, Gie, CHK03] for 
similar results for one dimensional random Schrodinger operators.) 

This allows us to prove linear Wegner estimates for random operators on very 
general metric graphs. In particular, we do not need any periodicity condition on 
the graph. Our assumptions are general enough to include, for instance, operators 
on the metric graph associated to I}" or, more generally, the Cayley graph of a 
finitely generated discrete group, or a Penrose tiling graph. 

Note that for quantum graphs there exists a certain dichotomy concerning the 
unique continuation principle. While along a single edge one has a quite strong 
version of this property, globally on the graph it does not hold. More precisely, a 
Laplacian (even without potential) on a metric graph with "free boundary condi- 
tions" at the vertices may exhibit compactly supported eigenfunctions. For ergodic 
Hamiltonians on discrete graphs it is well known (see for instance [DLM+03, KLS03, 
Ves05, Ves06, LV]) that compactly supported eigenfunctions are related to disconti- 
nuities of the IDS. Thus, for a complete understanding of the continuity properties 
of the IDS a systematic analysis of compactly supported eigenfunctions is necessary. 
See also [GLV] for related results and discussion. 

The structure of the paper is as follows: in the next section we formulate our 
results, Section 3 contains the bulk of the proof of the Wegner estimate, while 
arguments relating to boundary condition perturbations are deferred to Section 
4, and those to the unique continuation principle to Section 5. The last section 
contains a proof of the construction of the IDS via an exhaustion procedure. 

2. Model and Results 

First we define the geometric structure of the metric graphs and the operators 
acting on the associated L^-Hilbert space. 

Definition 1. Let V and E be countable sets and Q a map 

5: ^ X y X [0,cx)), e (i(e),r(e),/e)- 

We call the triple G ~ (y, E', a metric graph, elements oiV ^ ^(G) vertices, 
elements oiE = E{G) edges, t(e) the initial vertex of e, r(e) the terminal vertex of 
e and le the length of e. Both i(e) and T(e) are called endvertices of e, or incident 
to e. The two endvertices of an edge are allowed to coincide. The number of edged 
incident to the vertex v is called the degree of v. We assume that the degree is 
finite for all vertices. The elements of V which have degree equal to one we call 
boundary vertices and denote their set by . The elements oi -.^ V \ are 
called interior vertices. 

We will consider the graph G — {V, E, Q) as a topological space, more precisely as 
an one dimensional CW-complex, which we will again denote by G. The 0-skeleton 
of G is y and the collection of its one dimensional cells is given by E. Each one- 
dimensional cell e € E is attached either to one or to two zero dimensional cells 
V ^ V, namely t(e) and T(e). This defines the topological structure of G. 

We identify each edge e with the open interval (0,^e), where the point cor- 
responds to the vertex (-(e) and le to T(e). Assume that there exist constants 
< L , Z+ < oo such that for ail e £ E 

< I- < le < 1+ < oo. 
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The identification of edges by intervals allows us to define in a natural way the 
length of a path between two points in the topological space G. Taking the infimum 
over the lengths of paths connecting two given points in G, one obtains a distance 
function d: G x G ^ [0, oo). Since we assumed that each vertex of G has bounded 
degree, the map d is indeed a metric, cf. for instance Section 2.2 in [Sch06]. Thus 
we have turned G into a metric space {G,d). 

For a finite subset A C -E we define the subgraph Ga by deleting all edges 
e E E\A and the arising isolated vertices. We denote the set of vertices of Ga by 
Va, the set of vertices w £ VX of degree one (in Ga) by Vf, and its complement 
Va \ by Again, elements of are called boundary vertices of Ga and 
elements of interior vertices of Ga- Similarly as above we may consider Ga as a 
sub-CW-complcx of G with an induced topology and metric. 

For any A C E the Hilbert spaces L2{Ga) have a natural direct sum repre- 
sentation L2{Ga) = ffieGAi2(0, Ze)- In particular for A = i? we have L2{G) = 
©e6_EL2(0,/e), and for A C A C -E we have £2(G^) ^ ©egA-^2(0, ^e) C i2(GA) = 

©e6AE2(0,/e)- 

For a function (f>: G C and an edge e £ E we denote by := its restriction 
to e (which is identified with (0, Zg)). We denote by G(G) the space of continuous, 
complex- valued functions on the metric space (G,d). Similarly, G(Ga) denotes the 
space of continuous, complex- valued functions on the metric sub-space (Ga, d). For 
each V gV, any edge e incident to v, and function / e W^'^{e) C G^(e) = G^(0, le) 
we define the derivatives 

(1) def{v) := defiO) lim -^^"^ ~ -^^"^ if w = i{e) 
and 

(2) dJiv) -.^ dJik) lim /(^'-■) - /(^-^ - ^) if „ = ^(e). 

e\0 e 

Note that, since /|e G W^'^(e), by the Sobolev imbedding theorem the function / 
is not only continuously differentiable on the open segment (0, le), but also that its 
derivative has well defined limits at both boundaries and le- 

For any A C E it will be convenient to use the following Sobolev space 

W^^^(A) ©egAVt^'''(e) C C\A) := (BeeAC^e) 
with the norm ||0||^2,2(a) := W^^^Ww^^HoM' 

eGA 

Note that this space is defined on the edge set only and does not see the graph 
structure of G. The operators which we consider will be defined on the space 

V{Aa) -.^{f eW^'\A)nC{GA) I 

(3) ^v&Vl: J2 9<^f(^)= E ^-/W' V« G :/(«) = 0}. 

e^E ,L{e)—v e^E ,T{e)—v 

For each A C E we define a linear operator 

-AA:P(AA)^i'(A) 

by the rule 

(-Aa/)(.) 

if x G G is contained in the edge e. This way the function — Aa/ is defined on 
the set E C G, whose complement y = G \ in the metric space G has Hausdorff 
measure zero. 
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The operator — Aa is selfadjoint on the domain ©(Aa), see for instance [Kuc04, 
KS99] or Section 3.3 in [SchOG]. At the boundary vertices it has clearly Dirichlet 
boundary conditions, while the type of boundary conditions it has at the inte- 
rior vertices is called "free boundary conditions" by some authors and "Kirchhoff 
boundary conditions" by others. 

Next we describe the potential energy part of the Hamiltonian. Since it is ran- 
dom, we need an appropriate probability space. 

Let cj_ < cij+ be real numbers, f2 a probability space, and P a probability measure 
on il. Denote by E the mathematical expectation on fl with respect to P. Let 
W : X G ~> [oj-,uj+] be a stochastic process which is jointly measurable in the 
variables uj £ fl and x G G. For a fixed uj G H. we denote by W{uj) : L2{G) L-ziG) 
the multiplication operator W{uj)f{x) = {W{uj,x)f){x), x € G. 

To derive our Theorems 2 and 3 below we will need more specific hypotheses on 
P and W{u}). These are formulated in what follows. The first assumption describes 
random potentials of alloy type for which we are able to prove a Wegner estimate. 

Assumption 1. Let c+ > c_ > 0, s > and Cg be real numbers. For each edge 
e £ E let /ie be a probability measure which is absolutely continuous with respect to 
Lebesgue measure. More precisely, for e G E let g^^ £ Lac[i-iJ-,'-^+], \\ge\\oo £ Cg, and 
dfie{t) = 9e{t)dt on the interval [lj_,uj^]. Let the probability space Q, be given by 
the cartesian product Xeg£;[a;_,Cij+], and P by the product ^^eElJ-e- 
For each e G E let G ioo(0, le) be a single site potential satisfying 

C-XS, <Ue< C+XlO,l,], 

where Se C [0,le] is an interval of length \Se\ > s. We imbed Loo{0,le) — ioo(e) 
in Loo{G) and thus consider as an element of the latter space. Let the random 
potential W{uj) have the form 

WiLu):L2iG)^L2iG), = ^ WeUe- 

For quantum graphs with a Z'^-structurc we will establish the existence of the 
integrated density of states in Theorem 3. The following hypothesis formulates 
precisely what type of Z'^-structure we need for this result. 

Assumption 2. The vertex set V consists of the points 1,^ C M". The set of edges 
E consists of segments e = [a;, y\ parallel to the coordinate axes in where x,y G V 
have Euclidean distance equal to one. The union V U E C R'^ inherits the metric 
structure of . This structure coincides with the one defined earlier and we denote 
the corresponding metric space again by {G,d). 

Assume that for each x Gil' there is a measure preserving map t^; : J7 — > such 
that {tx}x£Z'' forms an additive group which acts ergodically on (f2,.7^, P). 

Assume that the random potential transforms in the following way under trans- 
lations: W{TxLU,y) = W{u!,y — x) for all x G 1^ G G and lo G Thus the 
Schrddinger operator is equivariant in the sense that 

(4) UxH{uj)U* = H{txUj) for all x G and uj G n 

where Uxf{y) = f{y — x),y G G denotes the unitary translation operator by x G 'L'^ . 

The restriction of the operator W{uj) to L2(A) will be denoted by Wa(^) or, if 
the set A is clear from the context, simply by W{uj) again. Note that the norm of 
the operator Wt^{ijj) is uniformly bounded in a; € and A C E. Thus the operator 
H\{lu) — Aa -I- Wa{uj) is selfadjoint on I>(Aa) for any ui G and A C E. If the 
set A C E contains a single element e we write — Ae(w) for — Aa(w) and He{oj) for 

A metric graph G together with a Schrodinger H operator which is defined on it 
we call a quantum graph. 
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Our main result is the following: 

Theorem 2. // Assumption 1 holds, then there exists for any A G M a constant C 
such that for all e G [0, 1] , 

IE{<^X[A-.A+e](i^AH)}<C-e-ttA. 

If in the terminology of Assumption 1 we have = [0, le\ for all e S -E, then the 
constant appearing in Theorem 2 can be chosen uniformly in the energy parameter 

A e M. 

Now we turn to the situation where the quantum graph has a Z^-structure as 
formulated in Assumption 2, and describe how the integrated density of states can 
be defined by an exhaustion procedure. 

For any I G N denote by A; the set of edges which are contained in 

{a; G I G (0, 1) for alH = 1, . . . , d}, 

and abbreviate Va, by Vi, and Ga, by G;. Denote by dGi the boundary of Gi 
as a subset of the metric space G, i.e. the vertices in Vi which have degree one. 
Similarly abbreviate A; := Aa, and Hiiuj) := H\^{uj). These are restrictions of the 
operators A and H{lu) to a finite cube with sidelength I and with Dirichlet boundary 
conditions on the boundary of the cube. The finite cube Schrodinger operator 
Hiiuj) is again self-adjoint, lower bounded and has purely discrete spectrum. Let 
us enumerate the eigenvalues of Hi{lo) in ascending order Xi{Hi{uj)) < X2{Hi{uj)) < 
X3{Hi{uj)) < . . . and counting multiplicities. 

Thus for each A G R and I G N, the counting function 

FtiX) := G N I A„(i/,H) < A} 

is monotone increasing and right-continuous, i.e. a distribution function, which is 
associated to a pure point measure. Denote by A'^(A) := jrF^{X) the volume-scaled 
version of F^{X). 

Theorem 3. Let Assumption 2 hold, then there exists a distribution function 
: M M and a subset fi' C o/ measure one such that for all uj £ and 
for all X where N is continuous the convergence 

(5) lim NUX) = NiX) 

holds. 

Remark 4. Under certain additional assumptions it is possible to prove that the 
normalized finite volume eigenvalue counting functions converge uniformly in the 
energy parameter to the IDS, see [GLV]. For metric graphs whose isomctry group 
does not exhibit a -structure, but is merely amenable, it should be still possible to 
define the IDS of an ergodic Hamiltonian by an exhaustion, proceeding analogously 
as in the paper [LPV04]. 

3. Proof of the Wegner estimate 

For the purposes of the proof it will be necessary to differentiate the spectral 
projection with respect to the energy parameter, which motivates the introduction 
of the following smooth 'switch function'. 

Let p be a smooth, non-decreasing function such that on (— cx), — e] it is identically 
equal to —1, on [e, cjo) it is identically equal to zero and ||p'||oo < Then 

X{x~.,x+e){x)<p{x~X + 2e)^p{x~X~2e)^ / p'{x-X + t)dt 

J-2e 
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Thus by the spectral theorem 



.2e 

X(A-.,A+e)(i?AM) < / dtp'{HA{u;)-X + t) 

in the sense of quadratic forms. Since -Be(A) = (A— e, A+e) is bounded and a{Hf^{Lo)) 
discrete, the above operators are trace class and we may estimate: 



tr 



XB,(\){Hf,{uj)) 



< tr 



2c 



2e 



riGN 



1 1'^" 



where A^(tj) denotes the eigenvalues of H/^iuj) enumerated in non-decreasing order 
and counting multiplicities. Only a finite number of terms in the sum are non-zero. 



Proof of Theorem 2. Let p be as above. Similarly as in [Kir96], p. 509, we estimate 

2e 

E{trx[A-c,A+c](i?AH)} < / ^ /'p'(A^H-A + Od<[]d/.e. 

To bound the right hand side we use Corollary 9 in Section 5 and the monotone 
convergence theorem to obtain the estimate 

[lo- ,C(J-|-J^^ —ze 

(6) < E / I[dpM ^ / l dp{Xt{.)-X+t) 



-2ca)_ 



Here we used the abbreviation A := A \ {e}. Denote by i?A(e,u'+) and H\{e,uj-) 
the operator 7?a, where the random variable Ue is set to its maximum respectively 
its minimum value. The n— th eigenvalues of this operators are abbreviated by 
A^(e,Lj_|-) and A^(e,a;-(_). Using monotonicity, the sum over n in (6) can be esti- 
mated by 

2c 

(7) llffellooE [ {pi>^n{e,^+)->^ + t)-piX^{e,LU_)-X + t)}dt. 

We introduce now the operators H'^{e,u!), * £ {D,N}, that coincide with Ha{ui) 
up to additional Dirichlet, respectively Neumann b.c. at the vertices i{e) and T(e). 
Their eigenvalues are denoted by A^'*. By Lemma 7 on Dirichlet-Neumann brack- 
eting in Section 4 we have 

(8) p(A^(e, LJ+)-X + t)- p{X^{e, c._) - A + t) 

< p(A^^(e, u;+)-X + t)- p(A^l'^(e, a;_) -X + t). 

Because of the decoupling of the edge e, the latter operators can be written as the 
direct sum H'^{e,uj) ~ H'^*(e,uj) H'^*{e,uj) of operators acting on £2(0, /e) and 
i2(GA\e), respectively. For the eigenvalues of these operators we use the notation 
A^'*(e,u;) and A^'*(e,(jj). Hence the sum over the terms in (8) can be separated in 
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the corresponding parts: 



(9) 



^ p(A^^(e, u;+)-X + t)- p(A^^(e,c._) - \ + t) 



(10) 



(11) 



+ E /'(^"''(e' - A + t) - p(A:^'^(e, - A + i). 



We estimate first the sum in (10). The difference iJ^ - ^^a 

is a perturbation of 

rank 2 in resolvent sense, see for instance [Sini95]. By Lemma 6 in Section 4 the 
first term in (10) obeys the bound 



where we used that the total variation of p equals one. Thus we are left with 
estimating 



Note that by definition of p we have p{x — X + t) < X[A-3e,oo)(a;) and ~p(x — X + t) < 
X(A+3£,oo) (a^) for all t S [— 2e, 2e]. Thus the trace in (12) is bounded by 

tr [X[A-3e,oo)(-ffA^(e,t^+)) - X(A+3e,oo)(-ffA^(e,t^-))] 

< tr[X[A-3£,oo)(-ffA^(e,W-) + - W_)||lie||oo) " X(A+3£,oo) (-^A ^ (e, ) )] 
= 2 + tr(x[A-3£-(a;+-c^_)||«e||^,A+3e)(^A^(e,t^-))) 

< C2. 

Here C2 is a constant which is independent of A and depends only on the considered 
energy interval. In fact, a careful look reveals that C2 depends only on the length 
(6e + (w+ — a;_)||we||oo) of the energy interval, but not on A £ M. Sec for instance 
Section 3 in [GLA']. 

Next we want to estimate the sum in (11). Let d := deg(t(e)) + deg(T(e)) — 2. 
Similarly as above, one sees that the difference i7^'^(e,cj+) — H'^^ {e,uj+) is a 
perturbation of rank d in resolvent sense. Consequently, the first term in (11) obeys 
the bound 



A telescoping argument bounds the whole sum in (11) by d times the total variation 
of p, i.e. by d. 

Hence the following upper bound on (6) completes the proof: 



p(A:;'^(e, C.+ ) - A + <) < p(A:f2(e,a;+) -X + t). 



By a telescoping argument we obtain the estimate 



E p( A^f 2 (e, a;+ ) - A + t ) - p( A^'^ (e, a;+) - A + < ) < 2 



(12) Pi^n'^ie, Lu+)-X + t)- p{X':^^{e, c._) - A + t) 



= tr [p(//l'^(e, u;+)-X + t)- p(Hl'^(e, c._) - A + t)] . 



p(A:;'^(e,c.+ ) - A + t) < p(A;;^^~(e,c.+) -X + t). 



Ci^lLgelloo / {d + C2)dt < 4CiiC2 + d)cg-e-iA^ 




□ 
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4. Changing boundary conditions 

In this section we show how to control the shifting of eigenvalues when one 
modifies boundary conditions. This is used in the proof of Theorem 2 above. 

Lemma 6 shows that changing the boundary conditions from the Dirichlet op- 
erator H^^ {e,uj^) to the Neumann one H^^ {e,uj-) shifts the eigenvalue index at 
most by the rank of the perturbation, i.e. by 2. The proof of this fact is based on 
the min-max formula for eigenvalues, see e.g. [RS78], Theorem XIII. 1. 

Lemma 7 is a monotonicity statement about boundary conditions known as 
Dirichlct-Neumann-bracketing. By quadratic form considerations one sees that the 
passage from Neumann to Dirichlet boundary condition shifts the eigenvalues up. 

First we proof an auxiliary lemma for bounded operators. 

Lemma 5. Let S,T ~ S+V be bounded, self adjoint operators on a separable Hilbert 
space TL with dim {ker V)^ = d < oo. Assume that the spectra of S and T below 
iuf aess{S), respectively mi u ess{T) , consist of an infinite, discrete set of eigenvalues. 

Then for the m-th eigenvalues of S and T, counted in ascending order including 
multiplicities, we have A„i(T) < Xm+diS) (Vm G N). 

Proof. Let L := (ker V)^. By the min-max principle one has 
A„i (T) = max min {(p \ Tip) 

dimL— m— 1 1, c/;_LL 

< max min _ [{ip \ Sp) + (p \ Vp))] 

dimL=m-l \\ip\\ = l. ip±L+L v ' 

=0 

< max min (p \ Sp) 

dhriL—m+d—l \\<p\\ — l, ipA^L 



□ 



The lemma above can of course not be applied to the Hamiltonians under consid- 
eration because they are unbounded operators. But there are some related operators 
for which the statement is true - one can simply compare suitable resolvents and 
rearrange the eigenvalues by the spectral theorem. 

Lemma 6. For the operators H'^^^ {e,uj-) and {e,uj^) we have 

(13) A™(i?l'^(e,c._)) < A„+2(i/l'^(e,c._)). 

Proof. We shift the spectrum of both operators by addition of a suitable constant 
and work in the following with the two arising strictly positive operators Hi (Dirich- 
let case) and H2 (Neumann case). 

Let Do := D{Hi) r\D{H2). Then one has 

(14) H^^ - \H2Da = Hi ^{H2 - Hi)H2^ |ff2-Do- 

By definition, i^ili),, ^ H2\do = 0, such that by (14) and continuity we get 

i.e. ker(Hf 1 - H^^)^ C H2D^ . 

We want to apply Lemma 5 to the bounded operators —H^^ and —H2^. So 
we have to show that dimker(iJ-f ^ — H^^)^ < 2, for what in turn dimiJ2£'(f < 2 
is sufficient. For this purpose we show next that Dq and D{H2) differ by a 2- 
dimensional subspace L. 

Let (f>i,(j>2 € C°°(0, ^e) be such that = 1 in a neighborhood of and 0i = in 
a neighborhood of 1^, and similarly let 02 = in a neighborhood of and (/)2 = 1 
in a neighborhood of le. Thus 0i, 02 are linearly independent vectors and moreover 
elements of D{H2) \ Dq. 
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Let -0 an arbitrary clement ol D{H2), ci ^/'(O) and C2 := ipHe)- Then 

(V' - C101 - C2(/)2)'(0) and (-0 - Cl(/)1 - C202)'(^e) = 

as well as 

(V' - Ci(/)i - C202)(O) = and - Ci^i - C202)('e) = 0. 

Hence ip — ci(j)i — C2(t>2 & Dq what implies that dim_D(_ff2) \ -Do = 2. 
Now, H2 : D{H2) -^H is one to one, and we get 

n = H2D{H2) 
= H2iDo + L) 
~ H2DQ + H2L. 

Hence dimiJ2£'^ < dimiJ2i < 2. 

So Lemma 5 is applicable to —H^^ and —H2^, and inequality (13) follows by 
the spectral theorem. □ 

It can be seen easily that the operators H'^^{e, w-) and iJ^'^(e, oj-) in the proof 
of Theorem 2 can be treated in the same way. 

Lemma 7. Let uj Cz ^l, A G E finite, and e S A 6e arbitrary. Consider the operators 
H\{e,LL!), H^{e,Lu) and Hf^{e,uj) and their eigenvalues X^{e,uj), A^'^(e,a;) and 
X!^'^{e,uj) (defined in Section 6). Then the following inegualities hold for all n £ N 

(15) A^^(e,c.)<A^(e,a;)<A^^(e,c.). 

Proof. Denote by h\{e,uj), h^{e,uj) and h^{e,uj) the quadratic forms associated 
to the operators HA{e,uj), H^{e,uj) and Hj^{e,uj). Then by Section 4.2 in [ScliOG] 
the quadratic form domains obey 

hf{e,uj) D /iA(e,a;) D h^{e,uj). 

Thus the statement of the lemma follows immediately if one applies the quadratic 
form version of the niin-max formula for eigenvalues, see e.g. [RSTs], Theorem 
XIIL2. □ 

5. Single site potentials of small support 

In this section we prove a uniform lower bound on the sum of derivatives of 
eigenvalues which is formulated in 

Lemma 8. 

(16) j:2|M>ft,/)>o 

eeA 

for all eigenvalues A^ of H\{uj) inside a bounded energy interval I . The bound C{I) 
does not depend on the set of edges A C E and on the eigenvalue index n G N. 

We infer immediately: 

Corollary 9. Let p : M ^ [0, 1] be a smooth, monotone function with p = —1 on 
{—00, — e] and p = on [e, 00). Then, for the n—th eigenvalue A^ of H\{ijj) we have 

p'(A:i(.)-A+t)<c,(/)EM^iM-A±^. 

eeA 

Proof. By the chain rule 

\ d^. = p(A„(a;)-A + 0)}_^^— , 

eeA eeA 
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(16) implies the estimate: 



p'{x^{co)-x + e)<c,ii)-'Y. 

eeA 



duJe 



□ 



To infer the lower bound (16), set S = UeeA and apply the Hellman-Feynman 
theorem, i.e. first order perturbation theory. For a normalized eigenfunction ipn 
corresponding to A^(a;) we have: 



E 

eGA 



dX^icu) 



duJe 



= y^ii'ri I Uelpn) > / \lpn\ 
eGA 



If the integral on the rightern side would extend over the whole of A, it would be 
equal to 1 due to the normalization of ^jj. A priori the integral over 5* could be 
much smaller, but the following Lemma shows that we can control the ratio of the 
two integrals. 

Lemma 10. Let I be a bounded interval and S,, C [0, 1^] a non-degenerate interval. 
There exists a constant Ci (/) > such that 



(17) 



\^V>c^{i) ivi 



for all A C E finite, all e G A and for any eigenfunction ip corresponding to an 
eigenvalue A G / of Ha{lu). 

Thus /g 1^/"^ > Ci{I) IV'P with the same constant as in (17). Hence Lemma 
10 implies directly Lemma 8. 

Proof. For y € M. denote by Se + y ■= {x G (0, /e) | x — y G Se} the translations of 
the set Se along the edge e. The derivative of the function 



(18) 



satisfies 



\ij{x)\^dx 



Sa+y 



\ip{x - y)\'^ dx 



■^ipix - y) 
dy 



d 



^{x - y) -^Tpix - y) dx 
dy 



iplx — y) dx - 

< '^U'\\LHS.+y)\H''\\LHS.+y)- 

Sobolev norm estimates (e.g. Theorems 7.25 and 7.27 in [GT83]) imply 

Wi^'h-'iS^+y) < Cs ML^S,+y) + H"\\L^S,+y)■ 

By the eigenvalue equation we have 



(19) 



(5.+,)=C4</>(y), Ci = C4\\W-X\U. 



Gronwall's Lemma implies (j){y) < exp(C4|y|) 0(0) and thus 



el JS. 



□ 
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6. Proof of Theorem 3 

The proof of Theorem 3 on the existence of the integrated density of states follows 
the arguments of [KM82]. There the convergence (5) of the IDS was established for 
all rational energies X € Q, whereas we establish the same fact for all A which are 
continuity points of the IDS. The main step consists in proving that a superadditive 
ergodic theorem of the paper [AK81] is applicable. 

In the whole of this section we assume that the quantum graph has a Z'^-structure 
as described in Assumption 2. 

Let us first describe the type of superadditive processes considered in [AK81]. 
Denote by T the semigroup of measurable transformations given hy t^, x £ Nq, 
where No := N U {0}, and by X the class of sets of the form 

{a; e R"^ I tti < Xi < hi, for all i = 1, . . . ,d} 

where a,b € Ng. For x e Nq and Q el denote by Q + x := {y e \ y - x £ Q} 
the translation of the set Q by x. 

Definition 11. A function F : I ^ L^{Vl,Y') which satisfies 

(i) Fq = Fq+x for all Q G X and a; e Ng, 

(ii) if Qi, . . . , Qn G X are disjoint sets and if their union Q ~ Uf^iQi is again an 
element of X, then 

n 

Fq>Y.Fq^, 

i=l 

(iii) 7(F):= sup ^E{Fq}<oo 

Qei,\Q\>o ' ' 

is called a (discrete) superadditive process. 

We state Theorem 2.4 from [AK81]: 

Theorem 12. If F is a discrete superadditive process and Qi :~ {x \ < Xi < I 
for all i = 1, . . . , d}, then 

lim^ 

l^oo l" 

exists for almost all cu £ il. 

In fact, in the case that the semigroup T acts ergodically on the probability space 
(ri,P) one can identify the limit, see the remark on page 59 in [AK81]: 

hm l^"^ FQ^{iL!) = 7(F) almost surely. 

/ — >oo 

In Section 6.2 of [KrcS.^)] the above statements are extended to the case that it is not 
the semigroup T which acts ergodically on Q, but rather the full group r^;, x G Z". 



To apply the superadditive ergodic theorem we consider for arbitrary, fixed A G R 
the process given by the eigenvalue counting functions of the Schrodinger operator 
with Dirichlet boundary conditions. For Q € X denote by A the set of edges 
e € E CM.'^ which are contained in the interior of Q and set 

Fq :=FQiX,Lo) :^i{n\Xn{HAiLu)) < A}, QeX. 

The Dirichlet Schrodinger operator Ha{lo) has been defined in Section 2. We will 
show that Fq, Q G X, is a superadditive process. 

Lemma 13. For any fixed energy value X, the process Fq is superadditive. 
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Proof. We have to check that the properties (i) ~ (iu) m Definition 11 hold. 

Since A = intQ n i? we have also A + a: = (intQ + x) n E for ah x e . 
The equivariance property (4) of the random operators carries over to the spectral 
projections and thus to the eigenvalue counting functions. Hence property (i) holds. 

Property (ii) can bee seen using Lemma 7 on Dirichlet-Neumann bracketing 
in Section 4. For Q and Qi,---,Qn £ I as in (ii), the set A contains A^ 
and finitely many edges ei, . . . , ejv which lie in int Q \ ljr=i Qi- introduc- 
ing Dirichlet boundary conditions at finitely many vertices one obtains from the 
operator Ha{uj) the direct sum 

n N 

(20) 0i/A,H © 0i^.,H. 

Hence the eigenvalue counting function of H/^{uj) is an upper bound of the one of 
the direct sum operator (20). Obviously the counting function of ®"^i Ha-{u}) can 
be estimated from above by the one of the operator (20). Now property (ii) follows. 

Denote the eigenvalue counting function of the negative Dirichlet Laplacian on 
e = (0, 1) by 

no{X) :^tt{"eN| A„(-Ae)<A}, 
Obviously the counting function of 

(21) -0Ae 

eGA 

equals |A| • no(A). Since the operators — Aa and (21) differ by Dirichlet boundary 
conditions at 2 |A|. or even less, vertices, the counting function of — Aa is bounded 

by 

|A| •no(A)+4d. |A|. 

For this conclusion we use an argument analogous to Lemma 6 in Section 4. 

Since the random potentials we are considering are uniformly bounded by a 
constant, say K, we have 

^q(A, w) < |A| • no(A + K) + Ad ■ |A| for ah A £ M, Q e I. 

The number of edges in the set A which is associated to a box Q e X is linearly 
bounded by the volume of Q. Hence (iii) is proven. 

□ 



Now we can complete the 

Proof of Theorem 3. For a fixed A S R one applies the ergodic theorem of [AK81] 
to the superadditive process Fq{X,uj), Q G X. Let us denote the corresponding 
j{F) by 7(A). By definition Fq{X,uj) < Fq{X,uj) for ah A < A G R and ah uj e n, 
Q G X, thus A I— > 7(A) is a non-decreasing function. In particular, it has at most a 
countable set of discontinuity points. Denote the complement of this set by C and 
choose a dense countable set S C C. Hence 7 is continuous at each A G S*. 

Since in our case the transformation group is ergodic, for each A there is a set 
fix of measure one on which the convergence lini;_>oo Fq^[uj) = 7(A) holds. 

Since S is countable, il' = Hags^^a still has full measure and the convergence 
statement of the superadditive theorem holds for all A G 5 and uj E ft' . Now define 
the distribution function A^(A) := lim^^^^^^ 7(A). Thus on the set C the functions 
7 and N coincide. 
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The monotonicity of A i-^ Fq^ (A, lo) and the continuity of N on the set C imply the 
convergence (5). To see this, choose a sequence A„ G S", A„ > A, hm„^oc An = 
Then we have 

r%,(A,c.) - N{X) < l-^FQ^Xn.io) ~ Af(A„) + iV(A„) - iV(A). 

For arbitrary w G fi' and e > we choose first n sufficiently large such that N{Xn) — 
N{\) < e/2 and then I sufficiently large such that l-'^FQ,{Xn,u}) - A^(A„) < e/2. 
Thus one sees that 

limsup 1~'^Fq,{X,u;) < N{X). 

I — *oo 

Similarly one can choose a sequence A„ G 5, A„ < A, lim„^tx) A„ = A and then show 
that liminf^^oo 1'^Fqi{^>^) > Thus the theorem is proven. □ 
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Note added: A few days before we finished this manuscript, the paper An- 
derson Localization for radial tree-like random quantum graphs by P. Hislop and 
O. Post, containing results not unrelated to ours, was made available at www.arxiv.org. 
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